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ABSTRACT 

On the basis of 16 years of spectroscopic observations of the four components of the gravitationally lensed broad absorption line 
(BAL) quasar H1413+117, covering the ultraviolet to visible rest-frame spectral range, we analyze the spectral differences observed 
in the P Cygni-type line profiles and have used the microlensing effect to derive new clues to the BAL profile formation. We first find 
that the absorption gradually decreases with time in all components and that this intrinsic variation is accompanied by a decrease in 
the intensity of the emission. We confirm that the spectral differences observed in component D can be attributed to a microlensing 
effect lasting at least a decade. We show that microlensing magnifies the continuum source in image D, leaving the emission line 
region essentially unaffected. We interpret the differences seen in the absorption profiles of component D as the result of an emission 
line superimposed onto a nearly black absorption profile. We also find that the continuum source and a part of the broad emission 
line region are likely de-magnified in component C, while components A and B are not affected by microlensing. Differential dust 
extinction is measured between the A and B lines of sight. We show that microlensing of the continuum source in component D has 
a chromatic dependence compatible with the thermal continuum emission of a standard Shakura-Sunyaev accretion disk. Using a 
simple decomposition method to separate the part of the line profiles affected by microlensing and coming from a compact region 
from the part unaffected by this effect and coming from a larger region, we disentangle the true absorption line profiles from the 
true emission line profiles. The extracted emission line profiles appear double-peaked, suggesting that the emission is occulted by a 
strong absorber, narrower in velocity than the full absorption profile, and emitting little by itself. We propose that the outflow around 
H1413-I-117 is constituted by a high-velocity polar flow and a denser, lower velocity disk seen nearly edge-on. Finally, we report on 
the first ground-based polarimetric measurements of the four components of H1413-H117. 

Key words. Gravitational lensing - Quasars: general - Quasars: absorption lines -Quasars: individual : H1413+ 1 17 

1 . Introduction stage, washing out their cocoons (Voit et al. 119931 Becker et al. 

I20TO . 

The broad absorption lines (BALs) observed in the spectra of ^^^^-^^ ^^^^ high-quality observational studies, in partic- 

quasars (or QSOs, quasi-stellar objects), blueshifted with re- ^j^^ ^^^^ spectropolarimetry (Ogle et al. [JMI, no clear view 

spect to the broad emission lines (BELs), reveal massive, high- ^j^^ geometry and kinematics of the BAL phenomenon has 

g . velocity outflows in active galactic nuclei (AGN). Such pow- ^^^^^^^ ^j^^j^ p^^^.^ spherically symmeti-ic winds appeared 

^ ; erful winds can strongly aff^ect the formation and evolution of ^^^^^^^ ^^e variety of observations (Hamann et 

^ . the host galaxy, enrich the intergalactic medium, and regulate ^^ ^jgg^ Qgj^ ^ggg^ equatorial disks, rotating winds, polar 

. ^ ■ theformation of the large-scalestructures (e.g. Silk and Rees ^^^^^ combinations thereof have been proposed more or less 

X] : Furlanetto ar^oeb mB Scannapieco and Oh |200l successfully to interpret the observations of individual objects 

Scannapieco et al.|2005]l. ^j. sm all gro ups of them (e.g. Murray et al [19951 Schm idt and 

_C3 ^ About 15% of optic aUy se lected quasars have BALs in their Hines 179991 Lamy and Hutsemekers l2004l Zhou et al. l2006l l. 

spectra (Reichard et al. 120031 ). Outflows may be present in all Given the large parameter space characterizing non spherically 

quasars if the wind is confined into a small soUd angle so that symmetric winds, BAL profile modeling must then be combined 

BALs are only observed when the flow appears along the fine with other techniques to determi ne the outflow properties in in- 

of sight (Weymann et al. I199U Gallagher et al. 12007b . On the dividual objects (e.g. Young et al llOOTl l. 

other hand, BAL QSOs could be quasars in an early evolutionary ^„ interesting method that can bring independent infor- 

mation on the quasars' internal regions is the use of gravita- 



* Based on observations made with the Canada-France-Hawaii ^lO"^! microlensing. Indeed, in a typical gravitationally lensed 

Telescope (Hawaii), with ESO Telescopes at the Paranal Observatory quasar, a solar mass star belonging to the lensing galaxy has 

(Chile) and with the NASA/ESA Hubble Space Telescope, and obtained an Einstein radius Re (the microlensing cross section) on the 

from the data archive at the Space Telescope Institute. ESO program ID: order of 10 ^ pc, which is comparable to the size of the con- 

074.A-0152, 075.B-0675, 081.A-0023. tinuum source. The microlens, moving across the quasar core 

** Mattre de Recherches du F.N.R.S. in projection, can successively magnify regions of area ^ 

Boursierdu F.N.R.S. inducing spectroscopic variations that could be used to extract 
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information on the quasar structure (Schneider et al. 119921 and 
references therein). Several studies, based on simulations, have 
demonstrated the interest of microlensing analyses for under- 
standing BAL QSOs ( Hutsem ekers et al . [19941 Lewis and Belle 
[19981 Belle and Lewis [20001 Chelouche [20051) . 

H1413+1 17 is a BAL QSO of redshift z ^ 2.55 showing typ- 
ical P Cygni-type profiles, i.e., profiles where the absorption is 
not detached from the emission. Turnshek et al. ( 1988 ) discussed 
the spectrum ofH1413+117in detail and made the first attempts 
to disentangle the emission from the absorption assuming an in- 
trinsic blue/red symmetry of the emission lines. H1413+117 is 
also a gravitationally lensed quasar constituted of four images 
(Magain et al. 119881 see Fig. [TJ. The lensing galaxy is faint 
and its redshift poorly known: indirect estimates give zi - 1-0 
(Kneib et al. [T99El l orz; ^ 1 .88 (Goicoecha and Shalvapin [20T0l l. 
Evidence of microlensing in component D has been suggested 
from both photometry and spectroscopy (Angonin et al. 11990! 
0stensen et al. |1997| l. In particular, Angonin et al. (1 19901 1 found 
that the equivalent width of the emission lines is systematically 
smaller in component D than observed in the other components, 
a result which can be interpreted by microlensing of the contin- 
uum source, the larger region at the origin of the emission lines 
being unaffected. This effect, which appeared to last at least a 
decade (Chae et al. 120011 Anguita et al. I2008I I. offers the possi- 
bility to separate the microlensed attenuated continuum (i.e. the 
absorption profile) from the true emission line profile, thus pro- 
viding new clues to the formation of BAL profiles (Hutsemekers 
[19931 Hutsemekers et al. [T994l l. 

In the present paper, we homogeneously analyze the spec- 
tra of the four components of H1413-I-117 obtained from 1989 
to 2005. The spectra cover the ultraviolet to visible rest-frame 
spectral range. In Sect. |3] we show that the spectral differences 
observed between the images can be consistently attributed to 
microlensing in spite of intrinsic variations. In Sect. |4| using a 
simple method, we separate the parts of the spectra affected and 
unaffected by microlensing, which basically correspond to the 
attenuated continuum and the emission lines. From these results, 
we derive a consistent view of the macro- and microlensing in 
H1413-I-117 (Sect. |5j- Finally, with the "pure" absorption and 
emission profiles in hand, we discuss the formation of the BAL 
profiles and the implications for the geometry and the kinematics 
of the outflow (Sect. O. 

2. Data collection 

Spectra of the four components of H1413+117 were gathered 
from archived and published data. Table [T] summarizes the char- 
acteritics of the spectra obtained over a period of 16 years, with 
the date of observation, the spectral range, the average resolving 
power R - /l/A/l, and the instrument used. 

The visible spectra secured in 1989 with the bidimensional 
spectrograph SILFID at the Canada-France-Hawaii Telescope 
(CFHT) are described in Angonin et al. ( I1990I I and Hutsemekers 
( |1993t . These spectra were obtained under optimal seeing con- 
ditions (0'.'6 FWHM). They provided the first spectroscopic evi- 
dence of microlensing in H1413+1 17. 

A series of spectra were obtained in 1993-1994 with the 
Hubble Space Telescope (HST) feeding the Faint Object spec- 
trograph (FOS). They cover the UV-visible spectral range (grat- 
ings G400H and G570H). These data are described in Monier 
et al. (I1998I I. A second series of HST spectra, yet unpublished, 
were obtained in 2000 using the Space Telescope Imaging 
Spectrograph (STIS) and the G430L grating (principal investi- 
gator : E. Monier; proposal # 8127). All HST data were retrieved 



Fig. 1. A deconvolved near-infrared image of the gravitationally 
lensed quasar H1413-I-117 with the four images and the lens- 
ing galaxy labelled (from Chantry and Magain [20071) . The im- 
age has been obtained in the F160W filter {A ^ 1.6 //m) with 
the NICMOS camera attached to the Hubble Space Telescope. 
North is up and East to the left. The angular separation between 
components A and D is 1" 1 . 



from the archive and reduced using standard procedures for long 
slit spectroscopy and prescriptions by Monier et al. ( I1998I I. 

In 2005, visible spectra were obtained with the integral 
field unit of the Visible MultiObject Spectrograph (VIMOS) at- 
tached to the European Southern Observatory (ESO) Very Large 
Telescope ( VLT). The data were obtained under medium quality 
seeing conditions (1'.'2 FWHM). Details on the observations and 
reductions are given in Anguita et al. ( I2008I I. For this data set, 
it was not possible to separate the spectra of images A and B of 
H1413+117. 

The 2005 near-infrared spectra obtained with the integral 
field spectrograph SINFONI at the VLT were retrieved from 
the ESO archive (principal investigator : A. Verma ; proposal 
075.B-0675(A)). Only the spectra obtained with the best see- 
ing (0'.'5 FWHM on May 22, 2005 in the H and K spectral 
bands) are considered here. The pixel size was 0'.'125x0'.'250 on 
the sky. The observations consist of four exposures per spec- 
tral band. The object was positionned at different locations on 
the detector for sky subtraction. The data were reduced using 
the SINFONI pipeline. Telluric absorptions were corrected us- 
ing standard star spectra normalized to a blackbody. The indi- 
vidual spectra were extracted by fitting a 4-gaussian function 
with fixed relative positions and identical widths to each im- 
age plane of the data cube using a modified MPFIT package 
(Markwardt 2009). Astrometric positions were taken from HST 
observations (Chantry and Magain [2007b . The spectra which ap- 
peared affected by detector defects and/or important cosmic ray 
hits were discarded. The good spectra were finally filtered to re- 
move remaining spikes. 

In addition, we have observed H1413+1 17 on May 10, 2008 
with the polarimetric mode of the Focal Reducer and low disper- 



D. Hutsemekers et al: Microlensing in the BAL QSO H1413+1 17 



3 



Table 1. Spectroscopic data 



Date 


Spectral range 


R 


Instrument 


(y/m/d) 








1989/03/07 


4400-6700 A 


~ 450 


CFHT + SILFID 


1993/06/23 


4600-6800 A 


1300 


HST + FOS 


1994/12/24 


3250-4800 A 


1300 


HST + FOS 


2000/04/21-26 


3000-5700 A 


500-1000 


HST + STIS 


2005/03/18 


3700-6700 A 


~ 250 


VLT + VIMOS 


2005/05/22 


1.95-2.45 


-4000 


VLT + SINFONI 


2005/05/22 


1.45-1.85 


~ 3000 


VLT + SINFONI 


2005/06/07 


1.10-1.40 /^m 


~ 2000 


VLT + SINFONI 



sion Spectrograph 1 (FORSl) installed at the Cassegrain focus 
of the VLT. Observations have been carried out with the Vhigh fil- 
ter, under excellent seeing conditions (0'.'6). Linear polarimetry 
has been performed by inserting in the parallel beam a Wollaston 
prism, which splits the incoming light rays into two orthogonally 
polarized beams, and a half-wave plate rotated to four position 
angles (e.g. Sluse et al. 12005b . In order to measure the polar- 
ization of the four images, the MCS deconvolution procedure 
devised by Magain et al. (11998b has been applied. We used a ver- 
sion of the algorithm that allows for a simultaneous fit of differ- 
ent individual frames obtained with the same observational setup 
(Burud l2001b . We constructed the PSF using a bright point-like 
object located ~ 15" fromH1413+l 17. Since this object is close 
to our target and is similar in brightness to the individual com- 
ponents of H14134-117, it provided a good estimate of the PSF. 
The Stokes parameters have been calculated from the photome- 
try of the quasar lensed images derived from the deconvolution 
process. 

3. Description of the spectra 

At the redshift ofH1413+117, spectra obtained in the UV- visible 
contain C iv ^^1548,1550, Si iv /1/11393,1402, N v ^^1238,1242, 
Pv ^^1117,1128,0 VI ^^1031, 1037,Lya ^1216 and Ly/3/ll026, 
while the near-infrared spectra contain Ha /16563, HyS /14861 and 
[Om] /t^4959,5007 (Figs. |2] and O. From the [Oiii] fines we 
measure the redshift z - 2.553. The UV resonance lines show 
typical P Cygni-type profiles with deep absorption while the 
Balmer lines show broad emission possibly topped with a nar- 
rower feature. Lyo' and Ly/? lines are weak due to absorption by 
the N V and O vi ions, respectively (e.g. Surdej and Hutsemekers 
fT987l l. 

3.1. Evidence for microlensing 

In Figs.|2]and[3] we compare the profiles of various spectral lines 
observed in the different images A, B, C, and D of H1413+1 17. 
A scaling factor is applied to superimpose at best the continua, 
considering in particular the continuum windows 4525^545 A 
and 5165-5220 A (i.e. 127 5-128 A and 1450-1470 A rest- 
frame, Kuraszkiewicz et al. 120021) . If the four images are only 
macrolensed, the line profiles observed in the spectra of the four 
components must be identical up to the scaling factor. On the 
other hand, line profile differences between some components 
may be indicative of microlensing, which is expected to mag- 
nify the -small- continuum region and not the -larger- broad 
emission line region. 

We first note that the line profiles in components A and B 
are essentially identical up to the scaling factor. This suggests 



that neither A nor B are strongly affected by microlensing. The 
small difference observed in the [O iii] lines may be related to 
the fact that the narrow line region could be partially resolved 
(Chantry and Magain l2007b . The wavelength dependence of the 
scaling factor reveals higher dust extinction along the B line of 
sight than along the A one, as discussed in detail in Sect. 13.21 In 
the following we consider the average spectrum of components 
A and B, denoted AB, as the reference spectrum unaffected by 
microlensing effects. As seen in the upper panel of Fig. [2] the 
spectrum of H1413+117 changes regularly with time, suggest- 
ing intrinsic variations in the quasar outflow (this is further dis- 
cussed in Sect. 13.31 ). These changes are observed in all compo- 
nents, indicating that the time scale of the intrinsic line profile 
variation is longer than the time delays between the four images. 
In fact, the longest time delay is not larger than one month ac- 
cording to the observations of Goicoecha and Shalyapin ( 12010b 
and in agreement with model^ (Kayser et al. 1 1 9901 Chae and 
Turnshek [T999] l. 

The spectrum of component D is clearly different when com- 
pared to AB. After scaling the continua, the emission appears 
less intense in D. This behavior is observed at all epochs and in 
the different spectral lines, superimposed onto the intrinsic time 
variations seen in all components. This is a clear signature of 
a long-term microlensing effect in D with amplification of the 
continuum with respect to the emission lines. The timescale of 
the effect is in agreement with previous estimates, i.e. on the or- 
der of 10 years (e.g. Hutsemekers |I993] l. As first pointed out 
by Angonin et al. (11990b . a difference is also observed in the 
absorption profiles. This difference is especially strong in the 
1989 and 1993 spectra. This is a priori not expected since the re- 
gion at the origin of the observed absorption lines has the same 
spatial extent as the continuum source. Differential microlens- 
ing of absorbing clouds smaller in projection than the contin- 
uum source has been proposed (Angonin et al. 11990!). However 
the timescale of such an event is expected to be smaller than 1 
year (Hutsemekers 1 1 993l l. ruling out this interpretation. Instead, 
we interpret the difference in the absorption profiles as due to the 
superposition of an emission line onto a nearly black absorption 
profile (SectlOTi. 

Although not as strong as in component D, spectral differ- 
ences are also observed when comparing C to AB. After scaling 
the continua, the emission lines appear slightly stronger in C at 
all epochs, with a possible differential effect between the broad 
Yip and the narrower [Oiii] lines. This suggests that microlens- 
ing also affects component C, de-amplifying the continuum with 
respect to the emission lines. 

The scaling factor used to fit the continuum of D to that one 
of AB clearly depends on wavelength. This can be explained ei- 
ther by chromatic microlensing (the source of UV continuum 
is less extended than the source of visible continuum and then 
more magnified), differential extinction (the extinction along the 
AB line of sight is higher than along the line of sight to D), 
wavelength-dependent dilution of the quasar continuum by the 
host galaxy light, or a combination thereof. Such a strong wave- 
length dependence of the scaling factor is not observed when 
comparing C to AB. The origin of this effect is further discussed 
in Sect. 15.21 



' We also estimate small time delays, typically less than 10 days, 
considering a classical singular isothermal ellipsoid (SIE)+shear model, 
a SIE+shear+galaxy model proposed by MacLeod et al. l l2009t . as well 
as using pixellated models with a symmetric mass distribution (Saha 
and Williams l2004] l. 
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Fig. 2. Intercomparison, at different epochs, of the Si iv and C iv line profiles illustrating the spectral differences between some 
images of H1413+117. Ordinates are relative fluxes. Vertical dotted lines indicate the positions of the spectral lines at the redshift 
z = 2.553 . The scaling factors needed to superimpose the continua are indicated. AB refers to the average spectrum of the A and 
B components (which could not be separated in the 2005 spectra). The upper left panel illustrates the time variation of the AB 
spectrum. 




Fig.3. Intercomparison of the Ly/? + Ovi and Lyor + Nv line profiles in the 4 images of H1413+117 at two epochs (lower panels). 
The comparison of the Hy6 + [Oiii] and Ha line profiles is shown in the upper panels. As in Fig.|2] ordinates are relative fluxes and 
vertical dotted lines indicate the redshifted line positions. The scaling factors needed to superimpose the continua are indicated. 
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< 




O.E 



12 3 

l/A(Mm) 

Fig. 4. The flux ratio Fa/Fb from the UV- visible spectra of 
1989, 1993, 1994, 2000, and from the near-infrared spectra of 
2005 (crosses of different colors). Photometric data point ob- 
tained at different epochs are superimposed (black diamonds 
with error bars). The continuous lines represent SMC-like ex- 
tinction curves redshifted to z/ = 1 .0. 



3.2. The extinction curve from A/B 

Since images A and B are not significantly aff'ected by mi- 
crolensing, the wavelength dependence of their flux ratio can 
be interpreted in terms of differential extinction in the lensing 
galaxy. Fig. |4] illustrates the observed flux ratio F^/Fb using 
all available spectroscopic data (slightly filtered and smoothed). 
Photometric data are superimposed. They were collected from 
Angonin et al. (119901) . 0stensen et al. (119971 die r atios F^/Fb 
are averaged per filter), Turnshek et al. (|1997| i. Chae and 
Turnsh ek (|20TO . Kneibet al. (fT998l l. Chantry and Magain 
(|2007] |. MacLeod et al. ( |2009] |. 

There is a lot of dispersion in the measured flux ratios which 
arises not only from inacurracies in the data but also from a 
possible combined effect of intrinsic photometric variations and 
time delay (cf. 0stensen et al. |1997l l. A clear trend is neverthe- 
less observed, indicating higher extinction along the line of sight 
to image B, as suggested by Turnshek et al. 1 1997 J . The flux ratio 
is reasonably well modeled using 



Fa 
Fb 



Fao 



'BO 



^Q-0.4fy)AAa 



(1) 



where AAb - Ab{A) - Ab(B) is the difference of extinction be- 
tween the A and B line of sights measured in the B filter and 
^{A) the extinction curve tabulated in Pei (119921) . redshifted to 
the lens redshift zi - 1.0 (Kneib et al. 119981 ). Since no obvious 
2200A feature typical of the Milky Way extinction is observed at 
2.3 yum ' (Fig. nil, we adopt a SMC-like extinction curve. A rea- 
sonably good fit is obtained with Fao/Fbo - 1.05 ±0.02 and AAb 
= -a09+0.02 (Fig.©. With z, ^ 1.88 (Goicoecha and Shalyapin 
12010] ) a similar fit is obtained with AAb = -0.055±0.015. For 
Zi - 1.0 (1.88), the lines of sight to images A and B sample re- 
gions of the lens galaxy located at respectively 4.9 (5.2) kpc and 
6.0 (6.4) kpc from the galaxy center, adopting the astrometry of 
Chantry and Magain (120071 ) and assuming a flat cosmology with 
Q„ = 0.27 and Hq = 70 km s ' Mpc '. The flux r atio Fao/Fbo 
= 1.19±0.10 measured by MacLeod et al. (12009b at 11 ;um in 
the mid-infrared, i.e. at wavelengths were both extinction and 
microlensing are expected to be negligible, is compatible with 



the extinction corrected flux ratio Fao/Fbo we derive, although 
marginally higher 

3.3. Intrinsic line profile variations 

Time variations in the absorption line profiles of BAL quasars 
are not uncommon (Barlow et al. 119891 119921 Gibson et al. 
120081 ). In H1413+117, Turnshek et al. (fT988l ) reported a deep- 
ening of the Si IV BAL between 1981 and 1985. On the contrary, 
between 1989 and 2005 (Fig. |2]), variations are observed as a 
gradual decrease of the depth of the BAL high-velocity part, the 
deepest component of the profile being essentially unaffected. 
Variations appear more complex in Si iv than in C iv, affecting a 
larger part of the absorption profile (see also Fig.|6|. 

The strongest change occurs between 1989 and 1993 and 
corresponds to an increase of the luminosity (Remy et al. 119961 
0stensen et al. 119971 ). Moreover, stronger absorption is accom- 
panied by stronger emission, which is an indication that reso- 
nance line scattering can play an important role in the emission 
line formation. In the Civ line of the AB spectrum, the high- 
velocity absorption appears ~15% larger in 1989 than in 2005 
while the emission is ~25% more intense. In the framework of 
resonance scattering where each absorbed photon is re-emitted, 
this may suggest that the high-velocity outflow has more scat- 
tering material perpendicular to the line of sight than absorbing 
material along the line of sight. This also requires a large cover- 
ing factor. 

4. Decomposition of the line profiles 

4. 1 . The method 



We follow a method similar to that used in Sluse et al. ( 120071 ). 
Assuming that the observed spectra F, of the different images 
are made of a superposition of a spectrum Fm which is only 
macrolensed and of a spectrum Fm^i which is both macro- and 
microlensed, it is possible to extract the components Fm and 
Fm^ by using pairs of observed spectra. Indeed, considering a 
line profile, we can write 

Fi = MFm + MfiFM^ (2) 
Fi = Fm + Fmh ■ (3) 

where M - M1/M2 is the macro-amplification ratio between 
images 1 and 2 and ju the micro-amplification factor of image 1 . 
We assume image 2 not affected by microlensing. To extract Fm 
and Fm^ when M is not known a priori, these equations can be 
conveniently rewritten 



Fm — 



-A 



Fmh - 



A-M 
M IF 



A-M 



(4) 
(5) 



where A - Mfi. Considering the line profiles as the sum of a 
continuum (absorbed or not in the blue) and an emission profile, 
we write F2 - F^ + F^ , Fi - F[. + F'^ , and then 



Fm — 



-A 
A-M 

Fc + 



F' F' 
A A 



M 



A-M 



F' 

— -F, 
M 



(6) 
(7) 



Up to a scaling factor, Fm only depends on A, while Fm^i only 
depends on M (Eqs. 4 and 5). A is the scaling factor between the 
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Table 2. Amplification factors determined from the (D,AB) pair 



Lines 


Date 




A 


M 






Si iv-C IV 


1989 


0.645 


±0.015 


0.390 


±0.020 


1.65 


±0.12 


Si iv-C IV 


1993 


0.915 


±0.015 


0.460 


±0.040 


1.99 


±0.21 


Si IV-C IV 


2000 


0.775 


±0.010 


0.400 


±0.020 


1.94 


±0.12 


Si IV-C IV 


2005 


0.775 


±0.010 


0.400 


±0.040 


1.94 


±0.22 


Lycf-N V 


1994 


1.050 


±0.025 


0.460 


±0.030 


2.28 


±0.20 


Lya-N v 


2000 


0.805 


±0.015 


0.430 


±0.030 


1.87 


±0.17 


U/3-[0m] 


2005 


0.630 


±0.010 


0.425 


±0.015 


1.48 


±0.08 


Ho- 


2005 


0.540 


±0.015 


0.395 


±0.015 


1.37 


±0.09 



Fi and F2 continua: A - F'^/Fc- It can be accurately determined 
as the value for which Fm{A) = in the continuum adjacent to 
the line profile, where F^ - F'^ - (Eq. 6). Assuming A con- 
stant across the line profile, Fm then only contains a contribution 
from the emission profile. Under the assumption that at least a 
portion of the observed emission profile is only macro-amplified 
(i.e. not micro-amplified), the macro-amplification factor can be 
estimated as the value of M such that FMfi{M) = F^ over that 
part of the line profile (Eq. 7, see also Appendix A). The micro- 
amplification factor of the continuum is then estimated at the 
wavelength of the line profile: = A /M. On the other hand, if 
the emission line is amplified exactly like the continuum, mi- 
crolensing cannot be distinguished from macrolensing. 

Eqs. 6 and 7 show that if the emission profile is only macro- 
amplified (i.e. if F'g - M F^), Fm/j only contains the underlying 
continuum Fc, while Fm contains the full emission profile F^. 
If the emission profile is micro-amplified to some extent, parts 
of it will appear in both Fmi, and Fm (see Appendix A). More 
specifically, (Fmii - F^) can be seen as the part of the emission 
profile which is not only macro-amplified, and Fm as the part of 
the emission profile which is not amplified like the continuum. 

The micro-amplification factor fj. and the macro- 
amplification factor M possess some specific chromatic 
behaviors. While // is assumed constant over the small wave- 
length range spanned by a line profile, it can be different 
at the wavelengths corresponding to different line profiles. 
Indeed, the micro-amplification of the continuum is related 
to the effective size of the continuum source which can be 
wavelength-dependent. The macro-amplification factor M 
may also contain a wavelength-dependent contribution due to 
differential extinction in the lensing galaxy, since extinction, 
like macrolensing, acts on the line profile as a whole. Finally, 
concerning the time dependence properties, M is expected to 
remain identical at different epochs of observation while fi can 
be time-dependent. 

4.2. The results 

4.2.1. Line profile decomposition from tine (D,AB) pair 

We extract Fm and Fm/^ from the spectra of H1413+117 using 
Fi = Fu and F2 = F\b in Eqs. 4 and 5. We have M < I and 
/i > 1 since D is fainter than AB and its continuum amplified 
(Sect. 13. 11 1. A is computed as the value for which Fm{A) = 
in the continuum windows adjacent to the line profiles (cf. 
Sect. 13. lb . Within the uncertainties, A is the inverse of the scal- 
ing factor determined in Figs. |2] and [3] M is computed as the 
largest value for which Fmii{M) > Fc over the whole line pro- 
file. For the resonance line profiles, the continuum may be com- 
pletely absorbed at some wavelengths such that Fm^{M) must 
also be larger than zero. Measurements of A and M are given 
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Fig. 5. The relation between A and measured from the spec- 
tral analysis of the (D,AB) pair (Table |2|. Diamonds represent 
data from Siiv-Civ lines, triangles data from Lytt+Nv lines, 
and squares data from the Balmer lines. The continuous lines 
represent A = MiJ where M - 0.39, 0.42, and 0.45. 



in Table |2] with jj. - A /M. Taking into account the noise and 
possible contaminating features in the observed spectra, a range 
of acceptable values is obtained which provides a rough esti- 
mate of the uncertainties. These are typically 2-3% for A and 
4-10% for M. As expected, the macro-amplification factor M is 
independent of the epoch of observation, the dispersion of the 
values being in agreement with the uncertainties. The temporal 
variation of A then essentially comes from the variation of the 
micro-amplification of the continuum. This is shown in Fig. |5] 
This figure also suggests that, within the uncertainties, there is 
no significant wavelength dependence of M due to a differential 
reddening between AB and D (see also Sect. 15. 11 1. 

With the measured A and M, we then compute Fm and Fm^, 
from Eqs. 4 and 5. The results of the line profile decomposi- 
tions are given in Figs. |6] to [8] Over large parts of the emis- 
sion profile (mainly in the wings), Fm/i - Fc. This indicates 
that at least a part of the observed emission line profile is unaf- 
fected by microlensing and that M actually represents the macro- 
amplification factor. We emphasize that the derived Fm and Fm^, 
profiles are robust against small changes of A or M. 

Fig. |6] illustrates the spectral decomposition for the Siiv - 
C IV region. Although different values of A (or fi) are found at 
different epochs, the extracted spectra are remarkably consistent. 
The microlensed part of the spectrum, Fm^, contains the contin- 
uum with the full absorption profiles as well as a small contri- 
bution from the core (not the wings) of the emission profiles. 
The bulk of the emission lines appears in the macrolensed-only 
part of the spectrum Fm, clearly showing a two-peak structure 
in C IV (at 5350 A and 5500 A). Recall that Fm shows the flux 
emitted from a large region of the quasar (much larger than the 
Einstein radius of the microlens) whereas Fm/i shows the flux 
emitted from a smaller region (comparable to and smaller than 
the Einstein radius). Interestingly enough, the small part of the 
emission profile observed in Fm/j is the core, i.e. the low-velocity 
part. The temporal variations of the absorption are particularly 
well seen in the Fm/^ spectrum of Si iv. 

The rest-frame UV spectra, and more particularly the Lya 
-t- N V region, are similarly analyzed (Fig. |7]i. Although the de- 
composition is less accurate due to structures in the continuum 
blueward of Lya (possibly due to narrow absorption features and 
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Fig. 6. The microlensed Fmii and macrolensed-only Fm spectra of H1413+117 extracted from the comparison of the D and AB 
spectra, at different epochs. The Siiv and Civ Hne profiles are illustrated. The amplification factors A and M are given for each 
epoch. Ordinates are relative fluxes, the Fmii spectrum being shifted upwards by 2 units. Vertical dotted lines indicate the positions 
of the spectral lines at z = 2.553. 
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Fig. 7. Same as Fig. |6] but for the Lye + N v line profile. 
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inaccuracies in the wavelength calibration), the extracted spectra Fm^i and Fm show the same quaUtative behavior as observed in 

the C IV and Si iv line profiles. 
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Fig. 8. Same as Fig.|6] but for the H/3 + [Oiii] (left) and Ho- (right) line profiles. 
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Fig. 9. Same as Fig.|6]and |2]but the microlensed Fm/^ and macrolensed-only Fm spectra are extracted from the comparison of the C 
and AB spectra obtained in 2000. The Si iv and C iv (left) and Lya + N v (right) line profiles are illustrated. 



Fig. [8] shows the decomposition for the H/? + [Om] and Ha 
line profiles. The micro-amplified spectrum is clearly a flat con- 
tinuum in the H/3 + [Om] spectral region while there is some 
evidence that the core (and not the wings) of the Ha emission 
line is micro-amplified. 

4.2.2. Line profile decomposition from tine (C,AB) pair 

The same kind of analysis can be done using the pair (C,AB). 
However the line profile differences are more subtle (cf. Figs.|2] 
and[3ll so that the value of M is closer to the value of A (i.e. fi 
closer to 1) in Eqs. 4 and 5, making the extracted spectra noisier. 
Only the spectra obtained in 2000, which show the most conspic- 
uous profile differences, are considered here. A de-magnification 
of the continuum explains the observations (Sect. 13.11 ) and the 
resulting Fm and Fm/^ are illustrated in Fig.|9l They are roughly 
similar to those derived from the pair (D,AB) but the part of 
the emission profile which is micro-amplified is different. While 
only a small part of the core of the emission is seen in the Fm/j 
profile derived from the (D,AB) pair, the red wing of the emis- 



sion profile is also observed in Fm^ computed from (C,AB). This 
red wing and the blue emission peak are not seen in Fm, suggest- 
ing that the high-velocity component of the C iv resonance line 
is micro-deamplified like the continuum. This is not unexpected 
since, for a given Einstein radius, demagnification regions with 
relatively smooth // variations can extend on larger scales than 
amplification regions (e.g. Lewis and Ibata I2004l i. The emis- 
sion line core, which appears in both Fm/j and Fm, should origi- 
nate, at least in part, from a region more extended than the high- 
velocity component. 

From the 2005 near-infrared spectra (the decomposition 
of which is not shown), we measure A = 0.74 for both the 
HyS+[Ora] and Ho' regions. The condition Fmij{M) > Fc is ver- 
ified at the wavelengths of Hp or Ha with M ^ 0.87, while 
FMfiiM) > Fe is obtained at the wavelength of the [Om] lines 
with M ^ 1.03. The latter value is comparable to the value 
of M derived from the UV-visible resonance lines. Although 
the narrow line region may be partially resolved (Chantry and 
Magain |2007| |. the [Om] emission lines are expected to origi- 
nate from a larger region and then less affected by microlensing. 
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making M ^ 1.03 a more plausible estimate. In this case, with 
// - A/M ^ 0.72 for the micro-amplification factor of the contin- 
uum, the Balmer emission lines do appear in both pM/t and Fm- 
This means that they are also micro-deamplified although not as 
much as the continuum. 

5. Lensing in H1413+117 

5.1. The macro-amplification factors 

In principle, the variation of M with the wavelength can be at- 
tributed to differential extinction. Unfortunately, for the (D,AB) 
pair, the wavelength dependence is not clear enough to extract 
an extinction curve, given the uncertainties on the determination 
of M (Fig.|5]and Table|2]i. The results nevertheless suggest that 
the differential extinction between AB and D is lower than be- 
tween A and B for which the extinction at Lya is ~ 1.2 times the 
extinction at Ha (Fig.|4]i. As a consequence, the value M(D,AB) 
= 0.395±0.015 determined at the wavelength of Ho', i.e. in the 
reddest part of our spectra, should not differ from the true macro- 
amplification factor by more than 2%. 

For the (C,AB) pair, we conservatively adopt M(C,AB) - 
0.95 ± 0.08 at the wavelength of Ha. Comparing with M(C,AB) 
= 1.15 ± 0.05 and M(C,AB) = 1.05 ± 0.05 at the wavelengths of 
Lyo' and C iv, respectively (Fig. |9]l, M(C,AB) might be slightly 
wavelength dependent, providing marginal evidence that extinc- 
tion is lower for C than for AB. Since the differential extinction 
remains low, we also assume that it does not affect the macro- 
amplification factor determined at the wavelength of Ha by more 
than 2%. 

The flux ratios with respect to component A are then 
Fbo/Fao = 0.95+0.02 (Fig. [Hi, Fco/Fao = 0.93+0.10 and 
Fdo/Fao = 0.39+0.04. The fact that Fco/Fao and F^o/Fao are 
different from the values estimated from photometry is due to 
a significant de-amplification of the C continuum and to a sig- 
nificant amplification of the D continuum, as derived from the 
analysis of the spectral lines. This emphasizes the need to prop- 
erly correct for microlensing before interpreting the flux ratios. 

MacLeod et al. ( |2009] l determined Fbo/Fao = 0.84+0.07, 
Fco/Fao = 0.72+0.07 and Fdo/Fao = 0.40+0.06 at 1 1 yL/m in the 
mid-infrared, i.e., where microlensing and extinction are thought 
to be negligible. Although only marginally different, the flux 
ratios of the B and C components relative to A seem slightly 
smaller than ours. If real, the origin of such a discrepancy is not 
clear but could be related to the intense starburst activity detected 
in the ho st galaxy of H1413+1 17 (Lutz et al. l2007l Bradford et 
al. I2009I I, which possibly contaminates with PAH emission the 
1 1.2 jjm (3.2 yum rest-frame) flux measurements. 

5.2. Microlensing of the continuum source 

The micro-amplification factor /i determined in component D 
depends on both the date and the wavelength (Table|2|i. 

Considering the Si iv-C iv spectral region, the strongest vari- 
ation of jj. occurs between 1989 and 1993 (see also Figs.|2]and 
|6]l. It roughly corresponds to a relative photometric variation be- 
tween A and D which can be observed in the V light curves of 
H1413+117 presented by Remy et al. d 19961 1 and 0stensen et 
al. (I1997I I, superimposed onto the common intrinsic variation of 
the 4 components. Between 1993 and 2000, the variation of jj is 
weaker, in agreement with the HST photometry in the F555W 
filter reported by Turnshek et al. (|1997| i and Chae et al. ( 120011 1. 

At a given epoch, yu decreases with increasing wavelength, 
suggesting chromatic magnification of the continuum source. 




_Q g I I I I I I I I I I I I I I I I I I I I I 

-1.0 -0.8 -0.6 -0.4 -0.2 -0.0 
log(Ar-est) 

Fig. 10. The micro-amplification factor // measured in 2005 
(Table IS against the wavelength of observation expressed in the 
quasar rest-frame (in fim). The straight line represents the model 
prediction (see text). 

This is best seen in the 2005 data (obtained within a 2 month in- 
terval) which span the largest wavelength range. We emphasize 
that /i, when determined from the line profiles, is not contam- 
inated by differential extinction (Sect. 14.11 ). In Fig. [TOl we plot 
the values of as a function of the wavelength of observation in 
the quasar rest-frame. 

The magnification /i of an extended source close to a caustic 
can be written 

M-l^o+ ^== m (8) 
yRs /Re 

where Rs is the source radius. Re the Einstein radius of the mi- 
crolens projected onto the source plane, ^ is a constant on the 
order of unity, ({d) a function which depends on the distance to 
the caustic and yUo a constant "background" magnification (e.g. 
Schneider et al. 1 19921 Witt et al. ll993b . In the framework a sim- 
ple model where the continuum is emitted by a Shakura-Sunyaev 
( 119731 ) thin accretion disk thermally radiating, Rs{X) oc A'^^^ (e.g. 
Poindexter et al. 12008b such that we may expect 

log (//-//(,) = -^log(^) + C. (9) 

As seen in Fig.[TOl this model nicely reproduces the data using C 
- -0.55 and assuming for simplicity /yq - 1 (i.e. no background 
(de- )magnification), supporting the idea of chromatic magnifi- 
cation of a continuum emitted by a Shakura-Sunyaev accretion 
disk. 

In principle, we could have used the Lya+N v and the C iii] 
/il909 emission lines, also present in the 2005 visible spectra, 
to measure n at other wavelengths. Although tentative estimates 
do agree with the observed trend, the quality of the data is not 
sufficient to derive reliable values of // at these wavelengths, due 
to the insufficient spectral resolution in the complex Lya+Nv 
region and to the fact that the Cm] line is truncated. Clearly, 
with better quality data, it could be possible to separate M and 
H at other wavelengths using additional line profiles and thus 
derive the temperature profile of the accretion disk. 

From the value of the constant C, we can derive a rough 
estimate of the size of the continuum source ; Rs{X)/Re < 
^4/3 IQ-2C jjjg Einstein radius is computed to be Re - 0.01 
\IM/Mq pc using zi-\ .0, a flat cosmology, Q.„, - 0.27 and Hq = 
70kms ' Mpc^'.Mis the mass of the microlens. Then, /?s(/l) < 
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Table 3. Polarimetry of the four images 





A 


B 


C 


D 


ABC 


p 


1999/03 


1.6±0.5 


2.3±0.5 


1.8±0.5 


2.9±0.5 


1.8±0.3 


1999/06 








1.0±0.6 


0.9±0.3 


2008/05 


1.4±0.1 


2.4±0.1 


1.2±0.1 


2.0±0.1 


1.65±0.04 


e 


1999/03 


75±9 


65+6 


71±8 


102±5 


70±4 


1999/06 








103±18 


87±11 


2008/05 


72±2 


79±1 


69±3 


96±2 


75±1 



The polarization degree p is given in percent and the polarization po- 
sition angle 9 in degree, East of North. The data obtained in 1999 are 
from Chae et al. ( 125011 . 



0.012 ^/M/Mq pc at the rest-frame UV wavelength /I = 0. 15 yum. 
This limits tightens to Rs(A) < 0.007 y/M/Me pc with zi = 1.88. 
These values agree with the radii obtained for the lensed quasars 
HEl 104-1805 and Q2237-H0305 on the basis of their photomet- 
ric variability due to microlensing (e.g. Poindexter et al. 120081 
Eigenbrod et al. l2008l Anguita et al. l2008bb . 

Dilution of the quasar continuum -microlensed- by the 
host galaxy light -not microlensed- can affect the interpre- 
tation of the wavelength dependence of fi. Denoting a - 
Fc(host)/Fe(qso) the ratio of the host and quasar continua at a 
given wavelength, we find that the micro-amplification factor of 
the quasar continuum /i^ is related to the measured fi by 

fig^fi(l+a)-a (10) 

if the host galaxy is unresolved and macrolensed as the quasar 
(i.e. if F'^ - Mfiq Fdqso) + MF dhost)), or by 

ti,^^iil+a)-a/M (11) 

if the host galaxy is resolved and not macrolensed (i.e. if F'^ - 
M fig Fciqso) + Fc(host) assuming the host contained in the mea- 
surement aperture). Dilution by the host galaxy is expected to 
be stronger in the rest-frame visible-infrared than in the ultravi- 
olet. In the rest-frame visible of a luminous high-redshift object, 
a < 0.15 (e.g. Schramm et al. 120081) so that fiq does not dif- 
fer from ju by more than 5% (even less if the host is partially 
resolved). This is negligible and within the uncertainties. 

5.3. Microlensing of a scattering region ? 

Chae et al. ( 1200 11 1 have obtained the first polarization mesure- 
ments of the four images of H1413-I-117 using the HST. The 
F555W filter was used. They noted that, in March 1999, the po- 
larization degree of component D might be higher that the po- 
larization degree of the other components (Table O. From this 
result, they suggested that microlensing also affects the scatter- 
ing region. The measurements obtained in June 1999 possibly 
indicate an intrinsic variation of the polarization observed in all 
components. Similar variations have been reported by Goodrich 
and Miller ([1995] ). 

Taking advantage of an excellent seeing, we were able to 
measure the polarization of the 4 components of H1413+117 
in the V filter, from the ground. Our measurements are also re- 
ported in Table [3] Within the uncertainties, the polarization de- 
gree of components A, B and C do agree with the March 1999 
values of Chae et al. (2001) , while the polarization degree of 
component D does not. Instead, we find that the difference be- 
tween the polarization degrees measured in A and C and those 



ones measured in B and D is significant. The difference between 
A and B is especially intriguing since we found no significant 
microlensing effect neither in A nor in B (at least before 2005). 
Possible interpretations could involve the polarization due to an 
extended scattering region resolved by the macrolens (possibly 
in the host galaxy, see Borguet et al. 12008b . or a differential po- 
larization induced by aligned dust grains in the lens galaxy. More 
data are clearly needed to correctly understand the meaning of 
these measurements. 

6. Consequences for the BAL formation 

In the previous sections we derived a consistent picture of mi- 
crolensing in H1413+1 17, showing that the continuum of com- 
ponent D (or more precisely all the regions of the quasar located 
in a cylinder of diameter ~ 2Re oriented along the line of sight 
and containing the continuum source) is magnified with respect 
to the more extended regions at the origin of the emission lines. 
This allowed us to disentangle the absorption part of the BAL 
profiles, essentially Fm^,, from the emission part, essentially Fm 
(Figs.|6]to[8]i. The observed profiles are equal to the sum of the 
Fm^i and Fm spectra (Eq. 3). The separation is robust against the 
uncertainties of the amplification factors. It is however not per- 
fect since emission which originates from regions close to the 
continuum source in projection appears in Fm/j- Selected spectra 
are illustrated in Fig.[TT]on a velocity scale. 

The absorption profile of the C rv BAL appears nearly black 
extending from v ^ -2000 km s"' to v ^ -10000 km s"'. It 
is especially interesting to note that the flow does not start at 
V = in the rest-frame defined by the [O m] emission lines. The 
part of the profile between -8000 and -10000 km s ' is clearly 
variable between 1993 and 2000, showing a smaller depth in 
2000. At a given epoch (1993-1994), the absorption appears 
stronger in N v and weaker in Si iv; this is best seen in the ve- 
locity range -2000 to -4000 km s ' and indicates an ionization 
dependence of the optical depth. The extracted emission profiles 
show a double-peaked structure which extends to the blue as far 
as the absorption profile does. The blue peak at -8000 km s"' 
appears much fainter than the red peak at -1000 km s The 
full emission profile (represented by the green line in Fig. [TTT i 
is roughly centered on the onset velocity of the flow (-2000 km 
s '), thus blueshifted with respect to the [Om] rest-frame. In 
fact, the full absorption + emission line profile appears in a rest- 
frame blueshifted by -2000 km s"' with respect to the rest-frame 
defined by the [Om] emission lines. Ha, on the other hand, ap- 
pears redshifted (Fig. [3]). Although not clearly understood, these 
line shifts are common in quasars (e.g. Corbin IT990I Mcintosh 
et al. 11999b and particularly strong in BAL QSOs (Richards et al. 
12002b . in agreement with our observations. Interestingly enough, 
Nestor et al. (12008b found a deficit of intrinsic Crv Narrow 
Absorption Line (NAL) systems at outflowing velocities lower 
than 2000 km s ', possibly due to overionization close to the 
accretion disk. 

The shape of the emission suggests that it is occulted by a 
strong absorber, narrower in velocity than the full absorption 
profile, and emitting little by itself. Very similar absorption and 
emission profiles are produced in the outflow model of Bjorkman 
et al. ( 11994b proposed for early-type stars. We build on this 
model to interpret our observations. A toy model, detailed in 
Appendix B, is used for illustrative purposes (a full radiative 
transfer modeling is beyond the scope of the present paper; it 
is presented in Borguet and Hutsemekers, 120101 where details 
on the flow geometry are also given). We assume that the out- 
flow in H1413+117 is constituted of two components: a quasi- 
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Fig. 11. Selected Fm (upper) and Fm^i (lower) spectra from Figs.|6]and|2l shown on a velocity scale. The zero velocity corresponds 
to the red line of the doublets redshifted using z - 2.553. Ordinates are relative fluxes. Fm is shifted upwards by 2 units. Vertical 
dotted lines indicate the positions of the spectral lines. The profiles from the illustrative model (Appendix B) are superimposed: 
the unabsorbed emission line (green), the disk absorption profile (red) and the disk-absorbed emission line (blue). Double-peaked 
emission profiles are produced using a disk absorption profile which is contained in the observed full (polar+equatorial) absorption 
line profile. 



spherically symmetric "polar" outflow, and a denser disk seen 
nearly edge-on. The equatorial disk expands slower than the po- 
lar wind and partly covers it. The polar outflow produces typ- 
ical P Cygni line profiles constituted of the superposition of a 
deep absorption extending from -2000 to roughly -10000 km 
s"' and a symmetric emission due to resonantly scattered pho- 
tons (e.g. Lamers and CassinelU |1999l ). This emission (assumed 
gaussian shaped for simplicity) is centered on v ^ -2000 km s"' 
and extends from -10000 to +6000 km s"'. Both the remaining 
continuum and the emission from the polar wind are absorbed in 
the equatorial disk. A double-peaked emission line is then pro- 
duced (Fig. [TTl i. little emission being expected from the edge- 
on disk. As we can see from Fig. [TT] this simple model is able 
to reproduce the main characteristics of the intrinsic emission 
line profiles extracted from the microlensing analysis. Since the 
disk is expected to also absorb the continuum, its absorption pro- 
file must be contained within the total polar+equatorial absorp- 
tion profile, as illustrated. Variability of the polar outflow optical 
depth will generate variations at the high velocity end of the ab- 
sorption accompanied by a change in the resonantly scattered 
emission, as observed (Sect. [331 see also Bjorkman et al.. 1 19941 
for simulations). Note that two-component winds have good the- 



oretical grounds (e.g. Murray et al. 119951 Proga and Kallman 
I2004I I and are supported by many observations interpreted with 
either disks or polar flows (cf. Sect.[TJ. 

Spectropolarimetric observations of H1413+117 (e.g. 
Goodrich and Miller 119951 Lamy and Hutsemekers 120041 ) 
provide additional evidence favoring this kind of scenario. 
First, the polarization angle rotates within the absorption line 
profiles, suggesting the existence of at least two sources and/or 
mechanisms of polarization. The polar outflow and the disk, 
expected to produce perpendicular polarizations, can play this 
role, especially in the case of BAL QSOs with P Cygni-type 
profiles (Goodrich 119971 Hutsemekers et al. 119981 Lamy and 
Hutsemekers I2004I I. Furthermore, the absorption in the polar- 
ized spectrum is clearly narrower than the absorption in the 
direct spectrum (this is best observed in Fig. 3 of Goodrich and 
Miller, |1995l l, supporting the existence of a slowly expanding 
equatorial disk which absorbs the polar-scattered flux. 

6. 1 . Microlensing in the BAL 

The difference observed in the BAL profiles of images AB and 
D (best seen in C iv, Fig.|2]i can also be interpreted in the frame- 
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work of this outflow model. In 1989 and 1993, the intrinsic ab- 
sorption at 5330 A (-9000 km s"') due to the high-velocity 
"polar" outflow is nearly black (Fig. |6] and [TTT i. In classical 
P Cygni line profile formation (e.g. Lamers and Cassinelli ll999l 
Hutsemekers and Surdei rT990t . this absorption is partially filled 
in with emission resonantly scattered at the same velocity (the 
blueward peak of the intrinsic emission, not absorbed by the 
slower disk, and seen in Figs. l6l and fTTTi. Since the emission line 
and the absorbed continuum react differently to the magnifica- 
tion by the microlens, a spectral difference is observed in the 
high-velocity part of the BALs seen in AB and D. Later, in 2000 
and 2005, the high-velocity part of the BAL profile is less opti- 
cally thick; the absorption is not as deep as in the nineties and 
the blue wing of the resonantly scattered emission which fills in 
the absorption is accordingly weaker (the blue emission peak ap- 
pears narrower or less intense in 2000 and 2005, Figs. l6l and [TTI). 
As a consequence, the microlens-induced spectral difference ob- 
served in the high-velocity part of the BALs appears smaller at 
these epochs. 



7. Conclusions 

Using 16 years of spectroscopic observations of the 4 compo- 
nents of the gravitationally lensed BAL quasar H1413+117, we 
derived the following results. 

- The strength of the BAL profiles gradually decreases with time 
in all components. This intrinsic variation is accompanied by a 
decrease of the intensity of the emission. 

- The spectral differences observed in component D can be at- 
tributed to a long-term microlensing effect, in agreement with 
previous studies. This effect consistently magnifies the contin- 
uum source of image D, leaving the broad emission line re- 
gion essentially unaffected. We also find that the continuum of 
component C is most likely de-magnified, while components A 
and B are not affected by microlensing. Differential extinction is 
found between A and B. 

- Using a simple decomposition method to separate the part of 
the line profiles affected by microlensing from the part unaf- 
fected by this effect, we were able to disentangle the intrinsic 
absorption (affected) from the emission line profile (unaffected). 
Consistent results are obtained for the different epochs of obser- 
vation. 

- Considering the macro- and micro-amplification factors esti- 
mated with this method, we obtain a coherent view of lensing in 
H1413H-117. In particular, we show that microlensing of the D 
continuum source has a chromatic dependence which is compat- 
ible with a continuum emitted by a standard Shakura-Sunyaev 
accretion disk. 

- To interpret the extracted absorption and emission line pro- 
files, we propose that the outflow from H1413+117 is consti- 
tuted of a high-velocity polar flow (at the origin of the intrin- 
sic variations) and a dense disk expanding at lower velocity and 
seen nearly edge-on. This is in agreement with spectropolarimet- 
ric data and supports the idea that BAL outflows can have large 
covering factors. 

In our analysis, we focused on the most robust results. 
Several interesting questions nevertheless remain opened, re- 
quiring more observations or a full radiative transfer modeling. 
In particular, why does the polarization of the different compo- 
nents apparently differ? Is the scattering region actually mag- 
nified? High-quality spectropolarimetry of the four images of 
H1413H-117 may solve this issue and bring more informations 
on the nature of the ouflow. As seen in Figs.|6]to|9] different parts 



of the emission line profiles appear micro-(de)amplified in com- 
ponents C and D. Detailed line profile calculations are needed to 
understand the origin of this micro-amplified emission and shed 
light on the outflow kinematics. With excellent signal-to-noise 
spectra, it would also be possible to better use the microlensing 
effect observed in component C. 
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Appendix A: Example of a line profile 
decomposition 

We consider a Hne profile constituted of an underlying contin- 
uum Fc(A) (absorbed or not at some wavelengths) and an emis- 
sion profile FgiA) = Ea{A) + Eh{A). We assume the continuum Fc 
micro-amplified by a constant factor yu^ , the component Eg of the 
emission micro-amplified by a constant factor /i^., and the com- 
ponent Eb unaffected by microlensing. We consider a typical 
case with yU^ > yU^. > 1 . If M is the relative macro-amplification 
factor between images 1 and 2, we have 



Fi = Mn,F, + Mn,E„+MEb 
f 2 = Fc + Eg + Eb ■ 



(A.l) 
(A.2) 



Using these expressions with A - Mfic and fi - jjc in Eqs. 4 
and 5, we find 



Fm — Eb + ^ En 



Mil 



Fc + ' 



Ea 



(A.3) 
(A.4) 



With He > He > Fm > and Fm^ > 0. As expected, Fm con- 
tains the emission profile Eb unaffected by microlensing, Fm/j 
contains the continuum, and both Fm and Fm^, contain a part of 
the micro-amplified emission profile E^. Up to a scaling factor, 
the micro-amplified profile £„ is given by Fm/i - Fc- 

To effectively compute Eqs. 4 and 5 and to determine the 
profile of FMfi, we need to know M0 (or fic)- Practically, we 
consider M as a free parameter in Eq. 5 and, by varying it, we 
adopt the value of M closest to A such that Fmh{M) > Fc at all 
wavelengths (Eq. 7). This is equivalent to adopt a value of p.c as 
close as possible to 1, thus ensuring that the macro- and micro- 
amplifications are best separated (see also Sluse et al. 120071 ). 



- It is always possible to write n'e^X) = E„ + Ej, where )ie = 
max(/j^(/t)). 

' As discussed in Sect. 14.11 the factor A and then the profile of 
are more easily determined. 



Unless M ^ A, this method provides a reasonably accurate esti- 
mate of M, and then of ^c- Indeed, denoting the free parameter 
M' and using expressions A.l and A.2 in Eq. 5, we write 



Muc -M' M-M' 

FMfiiM') - Fc = , E„ + — Eb 

^ A- M A- M 

such that FMfjiM') > Fc is satisfied when 

M' = MyUmin 



with 



min 



HeEg + Eh 

Eg + Eh 



over the line profile. 



(A.5) 



(A.6) 



(A.7) 



If Eg - and Eh ^ at some wavelengths, then fi^in - 1 and 
M - M'. This means that M can be derived empirically when 
at least a portion of the observed emission line profile Fc is not 
micro-amplified. In fact this is formally always true. Indeed, if 
Mmin + 1, the whole profile would be micro-amplified by [imm 
i.e. macro-amplified by Myimm instead of M, which is inconsis- 
tent with an optimal separation of macro- and micro-lensin^ It 
should nevertheless be kept in mind that if the emission line is 
amplified just like the continuum, microlensing cannot be distin- 
guished from macrolensing. Fortunately, it is realistic to assume 
that a portion -even a very small one- of the observed emis- 
sion line profile Fc is not micro-amplified given the large size 
of the emission line region compared to the source of the con- 
tinuum and to a typical Einstein radius. This hypothesis can be 
verified as follows: if we notice that Fm^i - Fc and Fm + 
over a certain wavelength range, then we must have M' - M 
and Hmm - 1- Indeed, combining Eqs. A.5 and A.6, we find that 
when FMfi - Fc then either Eg must be proportional to Eb over 
that wavelength range -which is unlikely since Eg and Eh are 
meant to originate from different emission regions- or yUmin - 1 ■ 



Appendix B: The absorption / emission toy model 

For a given line profile, we adopt for the disk optical depth 
the functional form 



T"d(w) = ^ To(/) exp ■ 



, w - w,.«(0 , 



(B.l) 



where w is the velocity. WccSS) is the position of the center of the 
absorption line /, and Wgg the width taken identical for all lines. 
The sum is computed over the 2 or 3 lines which constitute a line 
profile, with the appropriate weights ro(/). Similarly, we adopt 
for the unabsorbed emission line profile 



fe{w) = ^/rf)(0 exp- 



W - Wce(0 



(B.2) 



The disk absorption profile is computed as Fg - exp(-Tj) 
and the emission line profile absorbed by the disk as Fca - 
fc exp(-T(/). The total absorption of the continuum must be 
equal to the absorption in the polar flow multiplied by the ab- 
sorption in the disk. Although not computed, it can be obtained 
by inverting the observed intrinsic absorption profile. 

For all line profiles shown in Fig.[TT] we use Wea(l) = -5000 
km s-i, w,g = 800 km s"', w„(l) = -2000 km s"', w,, = 5000 



The amplification factors can be renormalized as follows: M' 
Mii^,^, yu;, = ^r/y"min and yu;'(/l) = yU^(^)/yUmi„. Notc that yu^i,, 
mm(ij!XX)). 
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km s ' , ; = 1 corresponding to the reddest line of the profiles, the 
position of the other ones being fixed by the doublet separation 
and/or by the Lya - N v velocity separation. to(/) and feo{i) are 
choosen to fit the observations. The parameters are not unique 
and may be different at other epochs. 



